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We present a measurement of pp — > — >■ {£ — e, fi) production with a data sample 

corresponding to an integrated luminosity of 6.2 fb~^ collected by the DO detector at the Fermilab 
Tevatron pp Collider. The results of the electron and muon channels are combined, and we measure 
the total production cross section and the differential cross section da/dp^, where p^ is the mo- 
mentum of the photon in the plane transverse to the beamline. The results obtained are consistent 
with the standard model predictions from next-to-leading order calculations. We use the transverse 
momentum spectrum of the photon to place limits on anomalous ZZ'y and couplings. 

PACS numbers: 12.60.Cn, 13.85.Rm, 13.85.Qk 



INTRODUCTION 

The standard model (SM) describes the electroweak 
interactions through a non-abelian gauge group SU{2)i^^ 
U{1)y, which includes self-interactions of gauge bosons. 
Because the Z boson carries no electric charge, a coupling 
between a Z boson and a photon is not permitted. The 
Z7 production in the SM is dominated by the lowest- 
order Feynman diagrams shown in Fig. [T] 
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An excess in the number of high-energy photons can be 
a sign of new physics, e.g., supersymmetry, as described 
in Ref. [T] or new heavy fermions with nonstandard cou- 
plings to the gauge bosons, as discussed in Ref. [2]. Such 
an excess of high-energy photons can be described by 
assuming only Lorentz and local U{l)em gauge invari- 
ant ZZ^ and Zjj trilinear gauge boson vertices of the 
form shown in Fig. [2] using an effective theory with eight 
complex coupling parameters, hY , where i = 1, 2, 3, 4 
and V = Z or 7 [3j. Here, the couplings parameters 
hY and hY {hY and hY) are associated with dimension- 
six (dimension-eight) operators which allow for an in- 
teraction between a Z boson and a photon. To con- 
serve tree-level unitarity at asymptotically high ener- 
gies, one can introduce form factors dependent on the 
square of the partonic center-of-mass energy, s, given by 
hY = /i^;/(l-l-s/A^)", where A is the mass scale at which 
the new physics responsible for anomalous couplings is 
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FIG. 1: Feynman diagrams for leading-order production 
in the SM: (a) and (b) initial-state radiation from one of the 
initial-state partons, (c) and (d) final-state radiation from one 
of the final-state leptons. 




FIG. 2: Feynman diagrams illustrating anomalous pro- 
duction with a ZZ'^ vertex (a) and a ^77 vertex (b) . 



introduced 31 ■ These anomalous gauge boson couplings 
would give rise to an excess of photons at high transverse 
momentum, pj, which can be searched for by measur- 
ing the total production cross section and the differential 
cross section da/dp7p for Z7 — >■ £+£"7 (££7 henceforth) 
production. If no evidence of new physics is seen, we can 
place limits on the real components of the CP-even cou- 
pling parameters, /103 and /104, for A — 1.2 and 1.5 TeV. 
Following Ref. i4j, we choose form- factor powers for the 
unitarity scaling dimensions of n = 3 for and n — A 
for h\ . Z7 production has been previously studied at 
collider experiments [5HT5]. and because the value of A 
greatly affects the scale of anomalous Z7 production, we 
choose to perform this analysis for the values of A that 
were used by the recent DO E] and CDF ^ analy- 
ses. This choice of A differs from the value used by the 
ALEPH [E', CMS [5], DELPHI [T(T, L3 TT], and OPAL 
fl2j collaborations. 

We present measurements of the inclusive cross sec- 
tion and differential cross section for Z7 production in 
the electron and muon channels using a data sample cor- 



responding to an integrated luminosity of 6.2 ± 0.4 fb~^ 
collected at -^s = 1.96 TeV by the DO detector at the 
Fermilab Tevatron Collider between June 2006 and July 
2010. These results provide a significant improvement 
in the sensitivity to anomalous ZZ^ and ^77 produc- 
tion compared to a previous DO publication, utilizing the 
same channels and an integrated luminosity of 1 fb~^ [S]. 
In addition to increasing the size of the data set, we also 
combine with a previous result in the same channels [S], 
along with another DO result [B] that used 3.6 fb~^ of 
Z7 — > vv^ production to place stringent limits on Z7 
anomalous couplings. 

THE DO DETECTOR 

The DO detector P^HT5] consists of a central track- 
ing system contained within a 2 T superconducting 
solenoidal magnet, surrounded by a central preshower 
(CPS) detector, a liquid-argon sampling calorimeter, and 
an outer muon system. The tracking system, consisting 
of a silicon microstrip tracker (SMT) and a scintillat- 
ing fiber tracker (CFT), provides coverage for charged 
particles in the pseudorapidity range |r/| < 3 [IH]. The 
CPS is located immediately before the inner layer of the 
calorimeter and has about one radiation length of ab- 
sorber followed by several layers of scintillating strips. 
The calorimeter consists of a central cryostat sector (CC) 
with coverage \r]\ < 1.1 and two end calorimeters (EC) 
which extend coverage to \ri\ « 4.2. The electromag- 
netic (EM) section of the calorimeter is segmented into 
four longitudinal layers (EMi, i = 1, 4) with transverse 
segmentation of Arj x A(/) = 0.1 x 0.1, except in EM3, 
where it is 0.05 x 0.05. The muon system resides beyond 
the calorimeter and consists of a layer of tracking detec- 
tors and scintillation trigger counters before a 1.8 T iron 
toroidal magnet, followed by two similar layers after the 
toroid. The coverage of the muon system corresponds to 
a pseudorapidity range \r]\ < 2. 

EVENT SELECTION 

Candidate Zj events are selected in the e~^e~j and 
fi'^fi'^j (ee7 and fxfi'y henceforth) final states. The pp 
interaction vertex must be reconstructed within ±60 cm 
of the center of the DO detector along the beam {z) axis. 
For the electron channel, a sample of candidate Z-boson 
events is collected with a suite of single-electron triggers. 
The electrons are selected by requiring an EM cluster 
in either the CC (It^I < 1.1) or EC (1.5 < \r]\ < 2.5) 
regions of the EM calorimeter with transverse momen- 
tum pt > 25 (15) GeV/c for the electron candidate with 
the highest (next-to-highest) transverse energy contained 
within a cone of radius A7^ = y/{Ar])^ + {Acj))^ = 0.2, 
centered on the axis of the EM shower. At least 90% 
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of the cluster energy must be deposited within the EM 
section of the calorimeter. Electron candidates, with a 
shower shape consistent with that of an electron, are re- 
quired to be spatially matched to a track and to be iso- 
lated in both the calorimeter and tracking detectors. To 
suppress jets and photons misidentified as electrons, a 
likelihood discriminant is built using a set of variables 
sensitive to differences in tracker activity and energy 
deposits in the calorimeter: the number of tracks and 
the scalar sum of the transverse momentum of all tracks 
within AT?, < 0.4 of the EM cluster, the fraction of energy 
deposited in the EM section of the calorimeter, the lon- 
gitudinal and transverse shower profile in the calorime- 
ter, and the ratio between the transverse energy in the 
calorimeter and the transverse momentum of the electron 
associated track. To further suppress jets misidentified as 
electrons, in particular, for high instantaneous luminos- 
ity conditions, a neural network algorithm is trained on 
DrcU-Yan Z/7* e'^e~ and jet data, using information 
from the calorimeter and CPS: the numbers of cells above 
a threshold in EMI within A7^ < 0.2 and 0.2 < AU < 
0.4 of the EM cluster, the number of CPS clusters within 
ATZ < 0.1 of the EM cluster, and the squared-energy 
weighted width of the energy deposit in the CPS. Events 
where both electrons are contained within the EC are ex- 
cluded because of the small signal acceptance. Candidate 
events where the Z boson decays into two muons are col- 
lected using a suite of single-muon triggers. Within the 
muon channel, muon candidates are required to be within 
I77I < 2 and matched to a well- isolated track in both the 
tracker and the calorimeter with transverse momentum 
Pt > 15 GeV/c. The highest px muon must have pt > 
20 GeV/c. Both muon candidates are required to origi- 
nate from within 2 cm of the interaction point in the z 
direction. 

Photon candidates in both the electron and muon 
channels are required to have transverse momentum p'^ > 
10 GeV/c within a cone of radius ATZ — 0.2 centered 
around the EM shower in the CC. The rapidity of the 
photon, T]^, is required to be \ri''\ < 1.1. Additionally, 
the photon candidate must satisfy the following require- 
ments: (i) at least 90% of the cluster energy is deposited 
in the EM calorimeter; (ii) the calorimeter isolation 
variable / = [£;tot(0.4) - £;EM(0.2)]/i;EM(0.2) < 0.15, 
where i?tot(0.4) is the total energy in a cone of radius 
ATZ — 0.4 and i?EM(0.2) is the EM energy in a cone 
of radius ATZ = 0.2; (iii) the energy-weighted cluster 
width in the EMS layer is consistent with that for an 
EM shower; (iv) the scalar sum of the pt of all tracks, 
p^™, originating from the interaction point in an annu- 
lus of 0.05 < ATZ < 0.4 around the cluster is less than 2.0 
GeV/c; (v) the EM cluster must not be spatially matched 
to either a reconstructed track or to energy depositions 
in the SMT or CFT detectors that are compatible with 
a trajectory of an electron [50]; and (vi) an output larger 
than 0.1 of an artificial neural network (Onn) [21^ that 



combines information from a set of variables sensitive to 
differences between photons and jets in the tracking de- 
tector, the calorimeter, and the CPS detector. 

The dilepton invariant mass, Mu, is required to be 
greater than 60 GeV/c^, and the photon must be sep- 
arated from each lepton by ATZ^^ > 0.7. Additionally, 
each lepton must be separated from a jet by ATZej > 0.5. 
In the electron and muon channels, we select 1002 and 
1000 data events, respectively. In order to reduce the 
contribution of final-state radiation (FSR), subset data 
samples are defined with the requirement that the recon- 
structed three-body invariant mass, Mu^, exceeds 110 
GeV/c^. With this additional requirement, 304 and 308 
data events are selected in the electron and muon chan- 
nels, respectively. 

Background subtraction 

The selected sample is contaminated by a small ad- 
mixture of Z+jet events in which a jet is misidentified 
as a photon. To estimate this background in the elec- 
tron channel, the fraction of jets that pass the photon 
selection criteria but fail either the Pj^"^ or the shower 
width requirement, as determined by using a dijet data 
sample, is parametrized as a function of p^ and rj^ (ratio 
method). The background from Z+jet production is then 
estimated starting from a data sample obtained by re- 
versing the requirements either on p^™ or on the shower 
width requirement, and applying the same parametriza- 
tion. A systematic uncertainty associated with the esti- 
mation of the number of real photons in the data sam- 
ples is due to the finite size of the dijet background sam- 
ple. After subtracting the estimated background from 
the data in the electron channel, we estimate 926 ± 53 
(stat.) ± 19 (syst.) signal events when no Mu^y require- 
ment is applied, and 255 ± 15 (stat.) ± 5 (syst.) signal 
events with Mu^ > 110 GeV/c^. 

To estimate the background in the muon channel, we 
use a matrix method to estimate the Z+iet background 
contribution. After applying all of the selection criteria 
described above, a tighter requirement on Onn is used 
to classify the data events into two categories, depending 
on whether the photon candidate passes (p) or fails (/) 
this requirement. The corresponding number of events 
compose a 2-component vector {Np, Nf). Thus, the sam- 
ple composition is obtained by resolving a linear system 
of equations (iVp,iV/)^ — £ x {Nz-y, Nzj)"^ , where Nz-y 
(Nzj) is the true number of Z -I- 7 (Z-|-jet) events in the 
fiducial region. The 2x2 efficiency matrix £ contains 
the photon Sj and jet Sjct efficiencies that are estimated 
using photon and jet Monte Carlo (MC) samples and 
validated in data. Based on these studies, the efficiencies 
are parametrized as a function of the photon candidates' 
77^ with 1.5% and 10% relative systematic uncertainties 
for and Sjot respectively. Having subtracted the esti- 
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mated background from data in the muon channel, we 
estimate 947 ± 40 (stat.) ± 16 (syst.) signal events when 
no Mu^ requirement is applied, and 285 ± 24 (stat.) ± 
2 (syst.) signal events requiring Mu^y > 110 GeV/c^. 

As a cross-check, the Z+jet background is also esti- 
mated through a fit to the shape of the Onn distribu- 
tion in data for both electron and muon channels, using 
MC templates constructed from simulated photon and 
jet events. The results are in good agreement with those 
obtained from the ratio and matrix methods. 



RESULTS 
Total cross section 

The total cross section for ££-f production is obtained 
from the ratio of the acceptance-corrected ££7 rate for 
Mu > 60 GeV/c2, AUi^ > 0.7, > 10 GeV/c, and 
1 77''' I < 1, to the total acceptance-corrected dilepton rate 
for Mu > 60 GeV/c^. Henceforth, these acceptance re- 
quirements are referred to as the generator-level require- 
ments. We utilize this method because uncertainties as- 
sociated with the trigger efficiencies, reconstruction ef- 
ficiencies, and integrated luminosity arc larger than the 
theoretical uncertainties and cancel in the ratio. This 
ratio is multiplied by a theoretical estimate for the to- 
tal cross section for inclusive Z/-f* — >• U production for 
M« > 60 GeV/c^: 



az-y X B 



. , ^31- X {gz X e)FEWZ- (1) 

Nff^- [A X e/D)« 



Here, Nf^^^ and Ng^^'^ are the number of measured Z and 
background-subtracted Zj events in the data sample, re- 
spectively. The factor {az x S)fewz calculated with 
the FEWZ next-to-next-to-leading-order (NNLO) genera- 
tor with the CTEQ6.6 parton distribution func- 
tions (PDF) [21]. The FEWZ theoretical prediction is 
262.9 ± 8.0 pb, where the dominant uncertainty is from 
the choice of PDF. The term B is the branching fraction 
for Z/Y U, which in the SM is 3.4% for either elec- 
trons or muons. The factor k corrects for the resolution 
effects that would cause events not to pass the selections 
on the generator-level quantities, e.g. a generator-level 
photon with < 10 GeV/c, but to pass the recon- 
struction requirements, e.g. a reconstructed photon with 
> 10 GeV/c. This factor is only used for — ^ ££7 
events, and corrects for the photon energy smearing that 
dominates in the first bin. The muon px resolution 
affects both Zf-y* ii and Z7 -J> U'y and the corre- 
sponding correction cancel in the ratio of cross sections. 
For the events that pass the generator-level requirements, 
the factors {A x €iD)u-y and {A x eiD)u provide the frac- 
tion of events that pass the analysis requirements, with 



all acceptances measured relative to the kinematic re- 
quirements at the generator level for the U and U'y fi- 
nal states, respectively. Events migrate between bins in 
p7p because of finite detector resolution, and these effects 
are taken into account in calculating [A x eiD)u'y as a 
function of p^^, while {A x e/£))« is calculated for the 
entire sample. To estimate k and A x e/£i, we use in- 
clusive Z/^* — >• U events generated with the PYTHIA |25] 
generator with final-state radiation simulated using PHO- 
TOS [22] and the CTEQ6.1L [57] PDF set. Because 
PYTHIA is a leading-order (LO) generator and does not 
reproduce the observed spectrum in data, generated 
events are weighted to reflect the distribution ob- 
served in Ref. [28]. Events are then traced through the 
DO detector using a simulation based on GEANT |29j . 
Data events from random beam crossings are overlaid 
on the simulated interactions to reproduce the effects 
of multiple pp interactions and detector noise. Simu- 
lated interactions that take into account the observed 
differences between data and simulation are reweighted, 
e.g., the z coordinate of the vertex, instantaneous lumi- 
nosity, trigger efficiency, lepton identification (ID) effi- 
ciency, photon ID efficiency, and resolution effects. Here, 
the factor {A x ej£i)u has values of 0.15 (0.17) in the 
electron channel (muon channel). When no constraints 
on Miij are applied, the factor k has average values of 
0.83 ±0.01 (stat.) and 0.85 ±0.01 (stat.) for the electron 
and muon channels, respectively, and the average value 
of {A X €10)11^ is 0.12 for both the electron and muon 
channels. Values for {A x e/£))«^ and k are similar for 
the subsample requiring Af«^ > 110 GeV/c^. 

To account for systematic uncertainty on the migra- 
tion into the sample from generated events with < 
10 GeV/c, we conservatively vary the number of events 
produced outside the generator-level requirements in the 
PYTHIA simulation by ±20%, found as an upper estimate 
in studies of photon energy resolution in this kinematic 
regime, to measure the effect on the final cross section 
measurement. We find that the effect introduces a 1.5% 
systematic uncertainty on the total cross section. The 
dominant uncertainty corresponding to the calculation 
of A X ejjj is due to choice of the PDF set. There are 
20 free parameters in the CTEQ6.1L parametrization of 
the PDF that reflect fits to data from previous experi- 
ments. The uncertainties on acceptance and efficiencies 
due to the PDF parametrization are estimated using the 
CTEQ6.1M PDF uncertainties, following Ref. |30j. We 
find a total PDF uncertainty of 3.5%, dominated by the 
uncertainty on the acceptance-correction to the full ge- 
ometrical lepton acceptance. The photon ID efficiency 
is determined from a simulated sample of photons and 
is estimated to have an uncertainty of 10% for p'^ < 15 
GeV/c and 3% for pj > 15 GeV/c. 

To reduce the contribution of FSR in the data sam- 
ples, we calculate the cross section with and without 
the M«-y > 110 GeV/c^ requirement. To combine the 
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electron and muon channels, we utilize the method in 
Ref. [31], which averages the results of measurements 
with correlated systematic uncertainties. We assume the 
PDF and photon ID efficiency uncertainties to be 100% 
correlated between the two channels. The total cross sec- 
tion results can be found in Tables U and |lll The measure- 
ments are consistent with the NLO mcfm [57 prediction 
using CTEQ6.6 PDF set |24j and the renormalization 
and factorization scales evaluated at the mass of the W 
boson, Mw = 80 GeV/cP. The PDF uncertainties as- 
sociated with the SM prediction are evaluated following 
Ref. innj. We reevaluate the values for the spectrum 
calculated by NLO mcfm with the renormalization and 
factorization scales set to 160 GeV/c^ and again at 40 
GeV/c^ and use these as estimates of the theoretical un- 
certainty of 1 standard deviation relative to the central 
NLO MCFM value. 



TABLE I: Summary of the total cross-section measurements, 
when no AIu^ requirement is applied, for individual channels, 
combined channels, and the NLO mcfm calculation with as- 
sociated PDF and scale uncertainties. 









az-y X B 


[fb] 


667 data 


1026 


± 


62 (stat.) 


± 60 (syst.) 


^fij data 


1158 


± 


53 (stat.) 


± 70 (syst.) 


combined data 


1089 


± 


40 (stat.) 


± 65 (syst.) 


NLO MCFM 


1096 


± 


34 (PDF) 


(scale) 



TABLE II: Summary of the total cross-section measurements, 
with the > 110 GeV/c^ requirement, for individual 

channels, combined channels, and the NLO mcfm calculation 
with associated PDF and scale uncertainties. 









az-y X B 


[fb] 


667 data 


281 


± 


17 (stat.) 


± 11 (syst.) 


/X/J7 data 


306 


± 


28 (stat.) 


± 11 (syst.) 


Uy combined data 


288 


± 


15 (stat.) 


± 11 (syst.) 


NLO MCFM 


294 


± 


10 (PDF) 


(scale) 



Differential cross section da/dplp 

We use the matrix inversion technique |33j to unfold 
the experimental resolution and extract d{az'y x B)/dp7p 
{da/dp7p henceforth), the differential cross section for 
Z"f — >■ U'^, as a function of the true p7p. The elements of 
the smearing matrix between true and reconstructed 
bins are estimated using the full simulation of the detec- 
tor response on a sample of Z7 events generated using 
PYTHIA. Then, the matrix is inverted to obtain the un- 
smeared spectrum. We confirm that the unfolding pro- 
cedure introduces a negligible bias. Following Ref. (34] . 



the position of the data points are plotted in Figs. [3] and 
[4] at the value of p7p where the cross section equals the 
average value for that bin. The theoretical uncertainties 
associated with the choice of PDF and the renormaliza- 
tion and factorization scales are determined analogously 
to the theoretical prediction for the total production cross 
section. The combined differential cross sections da/dp^ 
are shown in Figs. [3] and [4] for no Muj requirement and 
Mii^ > 110 GeV/c , respectively. The values associated 
with Figs. |3] and |4] are given in Tables |III| and |IV[ 
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FIG. 3: Unfolded da/dp^ distribution with no Mee-y require- 
ment for combined electron and muon data compared to the 
NLO MCFM prediction. 



LIMITS ON ANOMALOUS COUPLINGS 

To set limits on anomalous trilinear gauge boson cou- 
plings, we generate Z7 events for different values of the 
anomalous couplings using the NLO Monte Carlo gen- 
erator of Ref. U . SM Drell-Yan production is included 
by reweighting the spectrum to mcfm for vanishing 
anomalous couplings. As shown in Fig. [5] anomalous Zy 
couplings would contribute to an excess of high-energy 
photons as compared to the SM prediction. We apply 
the following generator-level requirements: Mu > 60 
GeV/c2, ATZi.^ > 0.7, p^ > 10 GeV/c, \ri"'\ < 1, and 
Mu^ >110 GeV/c^, generate p^ templates as a func- 
tion of the anomalous couplings, and use the known 
acceptance and resolution functions to fold the pre- 
dicted generator-level distribution into a reconstruction- 
level distribution for pj. Using Poisson statistics for 
p^ > 30 GeV/c, we define a likelihood function to com- 
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TABLE III: Summary of the unfolded differential cross section da/dp^, when no Mu^ requirement is applied, and NLO mcfm 
predictions with PDF and scale uncertainties 

M'y combined data NLO mcfm 

da/dp^ 
[fb/(GeV/c)] 

111.14 ± 4.40 (stat.) ± 11.99 (syst.) 104.02 ± 4.10 (PDF) tl.t (scale) 

51.41 ± 3.83 (stat.) ± 2.65 (syst.) 57.13 ± 2.23 (PDF) tH (scale) 

25.34 ± 2.74 (stat.) ± 1.13 (syst.) 28.77 ± 0.43 (PDF) tJ J, (scale) 

8.08 ± 1.45 (stat.) ± 0.40 (syst.) 10.16 ± 0.26 (PDF) tol (scale) 

3.23 ± 0.60 (stat.) ± 0.17 (syst.) 4.15 ± 0.16 (PDF) toil (scale) 

1.70 ± 0.26 (stat.) ± 0.088 (syst.) 1.60 ± 0.061 (PDF) (scale) 

0.34 ± 0.079 (stat.) ± 0.018 (syst.) 0.42 ± 0.017 (PDF) +°o°oll (scale) 

0.038 ± 0.014 (stat.) ± 0.002 (syst.) 0.052 ± 0.001 (PDF) toool (scale) 



bin 


center 


[GeV/c] 


[GcV/c] 


10-15 


12.4 


15-20 


17.2 


20-25 


22.5 


25-30 


27.5 


30-40 


34.4 


40-60 


48.5 


60-100 


76.5 


100-200 


124.5 



TABLE IV: Summary of the unfolded differential cross section da/dp^, with the Mee-, > 110 GeV/c^ requirement, and NLO 
MCFM predictions with PDF and scale uncertainties . 



combined data 



NLO MCFM 



Py bin 


pj center 


[GeV/c] 


[GeV/c] 


10-15 


13.7 


15-20 


17.2 


20-25 


22.0 


25-30 


27.4 


30-40 


34.5 


40-60 


48.6 


60-100 


76.5 


100-200 


124.5 



dcj /dpTp 
[fb/(GeV/c) 



13.48 ± 0.48 (PDF) (scale) 
12.25 ± 0.47 (PDF) t°oll (scale) 



13.57 ± 1.87 (stat. 
14.87 ± 2.17 (stat. 
7.91 ± 1.76 (stat.) 
5.30 ± 1.15 (stat.) 
3.08 ± 0.57 (stat.) 
1.73 ± 0.26 (stat.) 
0.34 ± 0.079 (stat. 



) ± 2.43 (syst.) 
) ± 2.30 (syst.) 
± 0.81 (syst.) 
± 0.44 (syst.) 
± 0.33 (syst.) 
± 0.17 (syst.) 
) ± 0.019 (syst.; 



0.038 ± 0.014 (stat.) ± 0.002 (syst. 



;.94 



± 0.25 (PDF) tofs, (scale) 
6.13 ± 0.21(PDF) to°2z (scale) 
3.71 ± 0.15 (PDF) 



-0.14 



(scale) 

1.57 ± 0.061 (PDF) + (scale) 
0.42 ± 0.017 (PDF) (scale) 
0.052 ± 0.001 (PDF) ttlll (scale) 
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FIG. 4: Unfolded da/dp'^ distribution with Mu-, > 110 
GeV/c^ for combined electron and muon data compared with 
the NLO MCFM prediction. 



TABLE V: Summary of the ID limits on the ZZj and Z77 
coupling parameters at the 95% C.L. 







U-y 6.2 fb" 
A = 1.2 TeV A 


1 

= 1.5 TeV 


Uj 7.2 fb"^ 
uwy 3.6 fb~^ 
A = 1.5 TeV 


"04 


< 
< 


0.050 
0.0033 




0.041 
0.0023 


0.026 
0.0013 




< 
< 


0.052 
0.0034 




0.044 
0.0023 


0.027 
0.0014 



pare the combined electron and muon channels with a 
predicted distribution for given values of anomalous cou- 
plings. In the absence of any significant deviation from 
the SM prediction, we set one-dimensional (ID) and two- 
dimensional (2D) limits on the anomalous coupling pa- 
rameter values at the 95% C.L. A combined log-likelihood 
function using all data is defined by the sum of the indi- 
vidual log-likelihood functions of the electron and muon 
channels. We include the effect of systematic uncertain- 
ties associated with transforming a Monte Carlo tem- 
plate from the generator-level into a reconstructed dis- 
tribution and find that these uncertainties contribute to 
the value of the calculated limits on the order of 1%. We 
generate a 10 x 10 grid of templates for the distri- 
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bution as a function of /ig^ and /104 for |/iq^| < 0.1 and 
< 0.01, while setting all other coupling parame- 
ters to zero, and the limits are derived by varying about 
the maxima of the log-likelihood functions for the 95% 
C.L [3S]. Results for the ID hmits for A 1.2 TeV 
and 1.5 TeV are shown in Table |Vl The ID and 2D hm- 
its on the anomalous coupling parameters are shown in 
Figs. [6] and [7] utilizing the electron and muon channels. 
In these figures, the dotted lines represent the theoretical 
limits on the anomalous coupling values, beyond which 
S'-matrix unitarity is violated. Because the /iq^ parame- 
ters come from dimension-eight operators, the limits are 
more constrained than those of Hq^ couplings, which are 
dimension-six. 



o 

C 10= 
n 

Q. 



D0, 6.2 fb ' • 



lly Data 
MCFM 

Scale uncert. 
PDF uncert. 
h;3 = 0-05,h; = 
A = 1.5 TeV 



Mi|^>110GeV/c^ 



10=^ 



p; [GeV/c] 



FIG. 5: The SM prediction and anomalous Z'j cou- 
pling production compared with the unfolded da/dp^ 
for combined muon and electron channels for > 30 
GeV/c and Muj > HO GeV/c^ 



We combine these results with those of a previous DO 
Z7 analysis [3] . In that analysis, the ID and 2D lim- 
its on the anomalous couplings parameter were calculated 
using a data sample corresponding to an integrated lumi- 
nosity of 1 fb^^ of data collected between October 2002 
and Feburary 2006 (3.6 fb~^ of data collected between 
October 2002 and September 2008) in the 667 and /i/i7 
channels (1/1/7 channel), for A — 1.5 TeV. Results can be 
found in Fig. |8] and Table |V| 



CONCLUSIONS 

We have measured the differential and total cross sec- 
tions for Z7 — ^ £ij production in pp collisions using the 
DO detector at the Tevatron Collider with and without a 
Mu^ > 110 GeV/c^ requirement. Both the total produc- 
tion cross sections and differential cross sections da /dpTp 
are consistent with the SM at NLO predicted by mcfm 
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FIG. 6: The 2D (contour) and ID (cross) limits on the anoma- 
lous coupling parameters for (a) ZZ'y and (b) ^77 vertices at 
the 95% C.L. for A = 1.2 TeV. Limits on S-matrix unitarity 
are represented by the dotted lines. 



|32) . We observe no deviation from SM predictions and 
place ID and 2D limits on the CP-even anomalous Z7 
couplings for A = 1.2 and 1.5 TeV. When combining 
with the previous DO analyses, the limits are compara- 
ble to those found in the most recent CDF result [7], 
which uses ~ 5 fb~^ in the and vv^ channels and 
A = 1.5 TeV. Our results include the first unfolded pho- 
ton differential cross section da/dp^, as well as the most 
precise measurement of the total production cross section 
of Z7 «7. 
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C.L. for A = 1.5 TeV. Limits on S-matrix unitarity are rep- 
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for A = 1.5 TeV. Limits on ^-matrix unitarity are represented 
by the dotted lines. 
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